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LONG-TERM  GOALS 

The  goal  is  to  investigate,  through  theory  and  by  analyzing  existing  data,  sea  surface  physics  and  air- 
sea  exchange  in  winds  that  range  from  weak  to  hurricane-strength.  Ultimately,  we  want  to  develop 
unified  parameterizations  for  the  fluxes  of  momentum,  sensible  and  latent  heat,  and  enthalpy  across  the 
air-sea  interface.  These  flux  parameterizations  will  provide  improved  model  coupling  between  the 
ocean  and  the  atmosphere  and,  in  essence,  set  the  lower  flux  boundary  conditions  on  atmospheric 
models  and  the  upper  flux  boundary  conditions  on  ocean  models. 

OBJECTIVES 

1 .  Develop  a  theoretical  framework  for  predicting  air-sea  fluxes  from  mean  meteorological  conditions 
and  apply  uniform  analyses,  based  on  this  framework,  to  datasets  that  we  will  assemble. 

2.  Assemble  a  large  collection  of  quality  air-sea  flux  data  that  represents  a  wide  variety  of  conditions. 

3.  Compute  fluxes  from  these  datasets  using  an  improved  analysis  that  better  accommodates 
measurements  made  over  heterogeneous  surfaces,  such  as  coastal  zones.  Focus  the  analyses  on 
common  problems  where  existing  bulk  formulations  perfonn  poorly — such  as  over  surface 
heterogeneity,  in  weak  winds,  and  in  very  strong  winds. 

4.  Develop  a  unified  algorithm  for  predicting  the  turbulent  air-sea  surface  fluxes  that  spans  the 
environmental  range  in  our  datasets,  obeys  theoretical  principles  and  constraints,  and  substantially 
exceeds  the  correlation  due  to  fictitious  correlation. 
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APPROACH 


Andreas  and  Mahrt  bring  diverse  expertise  to  this  project.  Mahrt  has  recently  been  focusing  on 
boundary-layer  processes  in  weak  winds,  when  stratification  and  surface  heterogeneity  are  important 
issues  and  when  Monin-Obukhov  similarity  theory  breaks  down.  Andreas,  in  contrast,  has  been 
concentrating  on  high  winds,  when  sea  spray  becomes  an  important  agent  for  modifying  the  usual 
interfacial  fluxes  of  momentum,  heat,  and  moisture.  In  additions,  Dean  Vickers  of  Oregon  State 
University  is  a  contractor  on  this  project.  Vickers  brings  expertise  in  processing  large  datasets — 
especially,  aircraft  data — and  in  parameterizing  air-sea  exchange.  Together,  we  will  develop  flux 
parameterizations  that  span  wind  speeds  from  near  zero  to  hurricane-strength. 


Traditionally,  the  air-sea  fluxes  of  momentum  (x,  also  called  the  surface  stress),  sensible  heat  (Hs), 
latent  heat  (HL),  and  enthalpy  (Qen)  are  fonnulated  as  down-gradient  fluxes  that  are  compatible  with 
Monin-Obukhov  similarity  theory: 


T  =  pu»  =  pCDrU^  , 

(la) 

H.  =  pcpCHrU,(0,[p-0,)  , 

(lb) 

Hl  =  pL.CE.Ur(Q,[t-Qr)  , 

(lc) 

Q,„  =  pCK[  Ur[cp(03fc  -  0r)  +  Lv(Qslfc  -  Q, )]  . 

(Id) 

Here,  p  is  the  air  density;  cp,  the  specific  heat  of  air  at  constant  pressure;  Lv,  the  latent  heat  of 
vaporization;  Ur,  the  wind  speed  at  reference  height  r;  0r  and  Qr,  the  potential  temperature  and  specific 
humidity  at  r;  and  0SfC  and  QsfC,  the  surface  temperature  and  specific  humidity.  Equation  (la)  also 
defines  the  friction  velocity  u*. 

The  crux  of  most  flux  algorithms  is  in  how  they  parameterize  the  transfer  coefficients  appropriate  at 
reference  height  r:  CDr,  CHr,  CEr,  and  Ckt  in  (1).  A  vast  amount  of  air-sea  interaction  literature 
describes  parameterizations  for  these  quantities.  In  this  project,  however,  we  take  a  different  approach. 
Andreas  (1992,  2010,  201 1;  Andreas  and  Emanuel  2001;  Andreas  and  DeCosmo  2002;  Andreas  et  al. 
2008)  has  suggested  repeatedly  that,  when  sea  spray  becomes  a  significant  agent  in  air-sea  exchange, 
(lb)-(ld),  in  particular,  are  not  accurate.  Moreover,  spray-mediated  exchange  becomes  “significant” 
(at  least  a  10%  effect)  at  modest  winds  speeds,  12-15  m/s. 

To  account  for  spray  effects,  we  formulate  the  total  scalar  fluxes  as 


Hs,T  -  Hs  +  Qs.sp  > 
HL,T  =  HL  +  Ql,sp  ’ 

Qe„,T  -  Hs  +  HL  +  Qen 


(2a) 

(2b) 

(2c) 


In  these,  subscript  T  denotes  the  total  flux  across  the  air-sea  interface.  The  first  tenn  on  the  right  [first 
two  terms  in  (2c)]  is  the  interfacial  flux,  parameterized  as  in  (1);  and  the  right-most  term  is  the  spray- 
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mediated  flux.  Because  the  spray-mediated  fluxes  do  not  scale  the  same  way  that  the  interfacial  fluxes 
do  [i.e.,  according  to  (1)],  the  traditional  transfer  coefficients  CHr,  CEr,  and  Ckt  are  inaccurate  for 
parameterizing  the  total  fluxes  Hs  T,  Hl,t,  and  Qen/r  in  high  winds  (Andreas  2011).  In  fact,  Andreas 
(2011)  recently  demonstrated  that,  when  spray-mediated  transfer  is  in  play,  sensible  heat  fluxes  can  at 
times  be  countergradient,  contrary  to  the  assumption  on  which  (lb)  is  based. 

To  address  these  various  ideas,  we  have  assembled — and  are  still  assembling — a  large  set  of  air-sea 
flux  measurements.  We  currently  have  in  hand  20  datasets  comprising  about  7000  air-sea  flux 
measurements.  In  this  set,  surface-level  wind  speeds  range  from  near  zero  to  72  m/s;  and  sea  surface 
temperatures  range  from  -1°  to  30°C.  This  dataset  thus  covers  almost  all  oceanic  conditions.  Flux 
parameterizations  developed  from  these  data  should,  indeed,  be  “unified.” 

WORK  COMPLETED 

We  investigated  the  variation  of  the  sea-surface  sensible  heat  flux  on  multiple  scales  using  data  from 
the  Gulf  of  Tehuantepec  Experiment  (GOTEX;  e.g.,  Romero  and  Melville  2010)  and  from  eight 
additional  aircraft  datasets.  We  applied  preliminary  quality  control  to  the  nine  datasets  and  discovered 
an  artificial  dependence  of  measured  wind  speed  on  aircraft  heading  in  several  of  the  datasets.  As  a 
result,  NCAR  has  reprocessed  the  GOTEX  data.  Using  these  aircraft  datasets,  we  have  also  carried  out 
a  study  of  the  variation  of  the  surface  heat  flux  and  drag  coefficient. 

Recently,  Foreman  and  Emeis  (2010)  suggested  that  the  drag  coefficient  formulated  as  (la)  is  ill-posed 
and,  instead,  tried  the  drag  relation 

u*  =  aUN10  +  b  (3) 


for  aerodynamically  rough  flow.  In  (3),  UNio  is  the  neutral-stability  wind  speed  at  a  height  of  10  m. 
On  the  basis  of,  perhaps,  a  thousand  data  points  from  the  literature,  Foreman  and  Emeis  confirmed  (3) 
and  found  a  =  0.05 1  and  b  =  -0. 14  when  both  u*  and  UNio  are  in  m/s. 

The  interesting  feature  of  (3)  is  that,  because  b  is  negative,  it  naturally  predicts  a  maximum  value  for 
the  neutral-stability,  10-m  drag  coefficient  (Cdnio): 


C 


DN10 


u* 


V  UNio  j 


a4JN10  +  b 


u 


N10 


=  a" 


1  + 


aU 


N10  J 


(4) 


Such  a  limited  drag  coefficient  in  high  winds  is  compatible  with  modern  hurricane  models  (Jarosz  et 
al.  2007;  Moon  et  al.  2007;  Sanford  et  al.  2007;  Chiang  et  al.  2011). 

Because  (3)  is  in  line  with  our  goal  of  developing  a  new  framework  for  parameterizing  air-sea 
interaction  that  relies  less  on  Monin-Obukhov  similarity  theory,  we  investigated  it  ourselves  with  much 
more  data  than  Foreman  and  Emeis  (2010)  had.  Figures  1  and  2  show  our  results. 

A  critical  part  of  developing  parameterizations  is  validating  them.  We  therefore  divided  our  datasets 
into  two  parts — primarily  chronologically  according  to  when  we  obtained  the  data.  Figure  1  shows  our 
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“original”  data,  which  come  from  ships,  towers,  and  aircraft;  Figure  2  shows  our  later  data,  which  are 
all  aircraft  data  obtained  generally  at  altitudes  less  than  40  m. 

In  both  Figures  1  and  2,  the  data  clouds  change  character  in  the  Unio  range  8-10  m/s.  For  wind  speeds 
below  this  range,  the  points  tend  to  turn  toward  the  aerodynamically  smooth  limit,  where  the  roughness 
length  is  given  as 


z 


Os 


v 

0.135—  . 
u* 


(5) 


Here,  v  is  the  kinematic  viscosity  of  air.  Hence,  in  both  Figures  1  and  2,  we  fitted  only  the  data  for 
which  UN10>9m/s  with  (3). 

In  Figures  1  and  2,  the  results  are,  respectively, 

u,  =  0.0581UN10  -0.214  ,  (6) 

u*  =  0.0584Un10  -  0.245  .  (7) 


Figure  1.  Eddy-covariance  measurements  of  the  friction  velocity  ufrom  our  “ original ”  dataset  of 

ship-board,  tower-based,  and  aircraft  data  are  plotted  against  the  neutral-stability,  10-m  wind  speed 
Umq.  The  red  line,  Eq.  (6),  is  the  best  fit  through  data  in  the  aerodynamically  rough  regime, 
UN10-9m/s.  The  green  line  shows  the  aerodynamically  smooth  limit.  The  gray  line  is  the  drag 

relation  that  Moon  et  al.  (2007)  inferred  from  their  hurricane  simulations.  The  orange  points  are 
the  aircraft  data  from  the  CBLAST  hurricane  flights  (French  et  al.  2007)  and  are 

ignored  in  our  fitting  of  the  red  line. 
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Figure  2.  As  in  Figure  1,  except  these  data  are  from  our  more  recent  set  and  were  all  collected  by 
low-flying  aircraft.  The  blue  line,  Eq.  (7),  is  the  fit  to  the  data  here  in  the  aero  dynamically  rough 
regime,  UN10  >9m/s,  and  is  almost  indistinguishable  from  the  red  line, 
the  fit  to  the  data  in  Figure  1 


In  both  of  these,  u*  and  UNio  are  in  m/s.  The  two  lines  are  almost  indistinguishable,  as  Figure  2  shows. 
In  other  words,  we  validate  (6)  with  (7),  or  vice  versa. 

RESULTS 

For  all  of  the  nine  “later”  aircraft  datasets,  upward  heat  flux  is  observed  for  slightly  stable  conditions 
(Figure  3).  The  magnitude  of  this  “countergradienf  ”  heat  flux  increases  with  wind  speed  and  is 
possibly  related  to  sea  spray  or  microscale  variations  of  surface  temperature  on  the  wave  scale. 
Mesoscale  heterogeneity  of  the  sea  surface  temperature  (SST)  can  also  produce  an  upward  area- 
averaged  sensible  heat  flux  for  slightly  stable  conditions  (Mahrt  and  Khelif  2010),  but  such  mesoscale 
variability  cannot  explain  the  general  increase  of  the  countergradient  heat  flux  with  increasing  wind 
speed. 

In  an  effort  to  reduce  offset  errors  and  different  SST  processing  and  calibration  procedures  among  field 
programs,  we  adjusted  the  SST  in  each  field  program  to  minimize  the  countergradient  flux  for  weak 
winds.  With  or  without  this  adjustment,  for  the  combined  dataset  in  Figure  3  (lower  panel),  the  extent 
of  the  upward  heat  flux  for  weakly  stable  conditions  increases  systematically  with  wind  speed. 
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Figure  3.  The  kinematic  surface  sensible  heat  flux,  w’6’,  from  4-km  flight  segments  for  all  nine 
aircraft  field  programs  in  our  “later”  dataset  is  composited  over  intervals  of  0sfc  -  0(z),  where  0sfc 

is  the  surface  temperature  and  €Kz)  is  the  potential  temperature  at  flight  level.  Upper  panel:  the 
basic  compositing,  where  the  10-m  wind  speeds  are  <  7  m/s  (green),  7-14  m/s  (black),  14-21  m/s 
(blue),  and  >  21  m/s  (red).  Lower  panel:  same  as  above  but  focusing  on  near-neutral  conditions. 
In  weakly  stable  conditions,  0sfc  -  0(z)<O°C,  the  sensible  heat  flux  is  countergradient  and 

increases  with  increasing  wind  speed. 
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In  neutral  stratification,  the  logarithmic  profile  gives  the  following  relationship  among  u*,  Unio,  and  Z(>: 


uN10  =  Yln(10/z«)  • 


(8) 


Here,  k  (=  0.40)  is  the  von  Karman  constant.  As  a  result,  we  can  estimate  the  roughness  length  (in 
meters)  in  aerodynamically  rough  flow  from  (3)  as 


z 


Or 


=  lOexp 


-k(u*  -  b) 
au, 


(9) 


In  turn,  we  can  predict  the  roughness  length  as  a  function  of  u*  for  any  aerodynamic  regime  by  adding 
(5)  and  (9): 


z0  =  0.135—  +  lOexp 
u» 


-k(u*  -  b) 


au, 


(10) 


This  result  also  yields  the  10-m,  neutral-stability  drag  coefficient  from  (8)  and  the  definition  in  (4): 


C 


DN10 


k 

ln(10/z0) 


|2 


(ID 


Figure  4  shows  (11),  where  we  have  used  the  coefficients  in  (7)  for  a  and  b. 
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Figure  4.  Various  expressions  for  the  neutral-stability,  10-m  drag  coefficient,  Como,  as  a  function 
of  Umo •  The  red  line  is  our  result,  (1 0).  The  blue  line  is  a  typical  Charnock  formulation.  The  gray 
line  is  the  Moon  et  al.  (2007)  result,  as  in  Figures  1  and  2.  Sanford  et  al.  (2007)  and  Chiang  et  al. 
(2011)  attempted  to  fit  the  results  that  Powell  et  al.  (2003)  deduced  from  dropsonde  profiles  in 
tropical  cyclones.  The  pink  symbols  are  from  Bell’s  (2010)  angular  momentum  budgets  for 
hurricanes  Fabian  and  Isabel.  His  results  are  so  scattered  that  they  do  not  help 
us  choose  among  the  various  models. 


Figure  4  also  shows  other  opinions  as  to  the  behavior  of  Cdnio-  The  “Charnock  +  Smooth”  line 
combines  the  Charnock  relation  for  zor  with  the  aerodynamically  smooth  limit,  (5),  and  is  common  in 
bulk  flux  algorithms  (e.g.,  Fairall  et  al.  1996).  Sanford  et  al.  (2007)  and  Chiang  et  al.  (2011)  attempted 
to  represent  the  Cdnio  values  from  Powell  et  al.  (2003)  in  their  hurricane  models.  Sanford  et  al.  used  a 
piecewise-continuous  fit  to  the  Powell  et  al.  data;  Chiang  et  al.  fitted  them  with  an  analytic  function. 

Moon  et  al.  (2007)  used  a  theoretical  model  for  wind-wave  coupling  and  flow  separation  in  their 
hurricane  simulations  and  inferred  the  roughness  length  from  their  modeled  surface  stress.  The  gray 
lines  in  Figures  1,  2,  and  4  represent  their  analysis.  They  never  presented  their  results  as  u*  versus 

UNio,  however,  and  therefore  probably  did  not  realize  that  their  theory  predicts,  essentially,  a  straight- 
line  relation  between  u*  and  Unio  for  Unio  above  about  20  m/s.  Mueller  and  Veron  (2009)  likewise 
produced  the  roll  off  in  Cdnio  by  modeling  just  wind-wave  coupling  and  flow  separation. 

Hence,  we  can  tentatively  conclude  that  wind-wave  coupling  and  increasing  wave  sheltering  is 
sufficient  to  explain  the  behavior  of  the  drag  coefficient  in  high  winds.  There  is  probably  no  need  to 
postulate  exotic  mechanisms  (mostly  involving  sea  spray)  to  explain  the  roll  off  in  Cdnio  with 
increasing  wind  speed  as  Barenblatt  et  al.  (2005),  Makin  (2005),  and  Soloviev  and  Lukas  (2010),  for 
example,  have.  Furthermore,  because  u*  seems  to  be  linearly  related  to  UNio  for  all  wind  speeds  in  the 
aerodynamically  rough  regime,  we  see  no  obvious  “drag  crisis”  (cf.  Ingel  2011). 

For  comparison,  Figures  1  and  2  show  the  CBLAST  aircraft  data  that  French  et  al.  (2007)  obtained  in 
hurricanes  Fabian  and  Isabel  in  2003.  In  all  cases,  the  aircraft  never  made  flux  measurements  lower 
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than  70  m,  and  half  the  runs  were  at  190  m  and  higher.  French  et  al.  tried  to  correct  the  flight-level 
momentum  fluxes  for  the  probable  flux  divergence  associated  with  measurements  this  high  but  seem  to 
have  been  not  entirely  successful:  In  Figures  1  and  2,  these  CBLAST  data  are  biased  low  compared  to 
the  measurements  nearer  the  surface.  We  have  therefore  not  included  these  data  in  fitting  (6)  and  (7). 
Bell  (2010)  estimated  Cdnio  from  the  angular  momentum  budget  of  an  axisymmetric  hurricane;  his 
data  came  from  dropsondes  released  in  hurricanes  Fabian  and  Isabel.  Clearly,  his  CDnio  values  are  no 
use  in  evaluating  which  of  the  drag  parameterizations  depicted  in  Figure  4  is  most  realistic. 

IMPACT/APPLICATIONS 

Our  inspection  of  the  aircraft  datasets  has  shown  widespread  problems:  The  measurements  depend  on 
aircraft  heading  with  respect  to  the  wind  vector.  We  are  continuing  to  work  on  this  problem.  Our 
analysis  also  indicates  that  current  fonnulations  of  air-sea  fluxes  based  on  Monin-Obukhov  similarity 
are  not  always  well  posed  for  data  analysis,  an  effect  that  can  mask  important  wave  effects. 

The  behavior  of  the  drag  coefficient  in  tropical  cyclones  has  been  a  crucial  knowledge  gap  at  least 
since  Emanuel  (1995)  reported  that  hurricane  models  could  not  produce  realistic  storms  if  their  drag 
parameterization  was  simply  an  extrapolation  of  drag  relations  obtained  at  lower  wind  speeds — like  the 
“Charnock”  relation  in  Figure  4.  Equations  (6)  and  (7)  now  provide  a  rational,  data-based  upper  limit 
for  Cdnio  in  hurricane-strength  winds.  Moreover,  our  analysis  explains  why  Cdnio  must  roll  off  with 
increasing  wind  speed:  There  is  no  “drag  crisis”;  rather,  there  was  a  crisis  in  understanding  how  to 
fonnulate  a  drag  law  in  high  winds.  Equation  (3)  seems  to  be  an  improved  approach. 

TRANSITIONS 

Journal  articles  and  conference  presentations  describe  our  work  on  air-sea  exchange.  Andreas  has  also 
developed  a  software  “kit”  that  contains  instructions  and  the  Fortran  programs  necessary  to  implement 
a  bulk  air-sea  flux  algorithm  that  includes  the  spray  effects  we  have  described  (see  Andreas  et  al.  2008; 
Andreas  2010).  The  current  version  of  that  code  is  3.4,  and  the  kit  is  posted  at 
http://www.nwra.com/resumes/andreas/software.php,  where  it  can  be  freely  downloaded. 

Andreas  is  currently  collaborating  with  hurricane  modelers  Isaac  Ginis  and  Tetsu  Hara  of  the 
University  of  Rhode  Island  to  get  these  new  parameterizations  for  the  air-sea  fluxes  into  their  hurricane 
model. 

RELATED  PROJECTS 

Andreas  is  in  the  second  year  of  a  three-year  project  funded  under  the  National  Ocean  Partnership 
Program.  This  project  is  on  “Advanced  Coupled  Atmosphere- Wave-Ocean  Modeling  for  Improving 
Tropical  Cyclone  Prediction  Models,”  with  Isaac  Ginis  at  the  University  of  Rhode  Island  and  Shuyi 
Chen  at  the  University  of  Miami  as  lead  Pis.  Andreas  is  a  subcontractor  to  the  University  of  Rhode 
Island  and  will  supply  expertise,  code,  and  analyses  to  help  the  project  understand  how  sea  spray 
affects  hurricane  intensity. 
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